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We have investigated the antiferromagnetic phase of the 2D, the 3D and the extended Hubbard 
models on a bipartite cubic lattice by means of the Composite Operator Method within a two-pole 
approximation. This approach yields a fully self-consistent treatment of the antiferromagnetic state 
that respects the symmetry properties of both the model and the algebra. The complete phase 
diagram, as regards the antiferromagnetic and the paramagnetic phases, has been drawn. We firstly 
reported, within a pole approximation, three kinds of transitions at half-filling: Mott-Hubbard, 
Mott-Heisenberg and Heisenberg. We have also found a metal-insulator transition, driven by doping, 
within the antiferromagnetic phase. This latter is restricted to a very small region near half filling and 
has, in contrast to what has been found by similar approaches, a finite critical Coulomb interaction 
as lower bound at half filling. Finally, it is worth noting that our antiferromagnetic gap has two 
independent components: one due to the antiferromagnetic correlations and another coming from 
the Mott-Hubbard mechanism. 
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I. INTRODUCTION 

Since almost half of a century the Hubbard model^ is 
one of the fundamental models in the theoretical descrip- 
tion of strongly correlated electron systems. Despite its 
simplicity it has always been taken as one of the proto- 
types when modeling electrons in narrow bands. Within 
this context it had wide applications in the theory of 
magnetism^^^ and metal-insulator transitions.^ The dis- 
covery of the high- Tc superconducting cuprate materi- 
als,^ where the special properties are to a large extent 
due to strong correlations,^ gave a new impetus on the 
investigation of the Hubbard model and its derivates as 
the t-t'-U, the extended Hubbard and the t-J models. 
To gain a better theoretical understanding of the high- Tc 
materials, not only the superconducting phase of the pro- 
posed models is studied, but also their normal phase and 
such with further broken symmetries. The antiferromag- 
netic state is of particular interest, as it is believed that 
the pairing mechanism in the cuprates is mainly due to 
magnetic correlations in an itinerant electronic system. 
In addition, the high- Tc compounds show a wide range 
of metal-insulator transitions; in particular, a phase tran- 
sition from an antiferromagnetic insulator to a supercon- 
ductor is observed. The antiferromagnetic phase of the 
Hubbard model has been subject of intensive study by 
both numerical and analytical methods. 

Despite the simplicity of the model there are no ex- 
act solutions known except for a few special cases^'^ and 
one must recur to approximate solutions. The Dynamical 



Mean Field Theory (DMFT)^ and a variety of projection 
techniques^°"^^ are among the most popular approaches 
to the Hubbard model and its derivates. 

The Composite Operator Method {COM)^^'^^'^^ that 
we present here belongs to the above mentioned class of 
projection techniques. Choosing as elements of the fun- 
damental spinor a set of field operators built up from the 
electronic ones (hence the name), the equation of mo- 
tion for the Green's function of this spinor are obtained. 
These equations have been solved by means of both a pole 
approach^^' and a two-site resolvent method. The 
use of composite operators asks for a careful choice of 
the Fock space where the Green's function is realized. 
We have shown that the Pauli principle plays a funda- 
mental role in unambiguously fixing the representation 
by determining the parameters that appear in the scheme 
owing to the non-canonical algebra of the composite op- 
erators. Let us note that by Pauli principle we mean 
all relations among operators dictated by the algebra. A 
comprehensive treatment of the paramagnetic phase of 
the Hubbard model within the two-pole approximation 
has led to a good agreement with both numerical studies 
and some experimental properties of the cuprates. In 
the case of the ID Hubbard model we have shown that 
the two-pole approximation reproduces almost exactly 
the Bethe ansatz results for the ground state energy. 

In the last years an intensive study by the Spectral 
Density Approach (SDA)^^'^"^ of the magnetic phases of 
the Hubbard model on different lattice types and by vary- 
ing dimensionality has been done^''"^^ and the complete 
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magnetic phase diagram of the Hubbard model on a 3D 
fee-lattice was derived. The SDA corresponds to an ex- 
pansion of the Green's function in terms of the electronic 
spectral moments. It has been shown that this approach 
corresponds to the COM when a specific choice of the ba- 
sic spinor is made and the pole approximation is used.^^ 
However, this correspondence is only relative to the func- 
tional dependence of the Green's function and huge dif- 
ferences arise when the representation is not properly 

fixed.20.21,28 

To complement the analysis of the paramagnetic phase 

for the Hubbard model in the two-pole approximation by 
the COM, wc here present a study of the antiferromag- 
netic phase within this framework. As the most simple 
case of antiferromagnetism, where in general the spin and 
the translational invariance arc broken, we consider an 
antiferromagnetic state on a bipartite lattice character- 
ized by its staggered magnetization. It turns out that 
in this phase the simple Hubbard model does not pro- 
vide enough internal parameters to fix the representation. 
We therefore investigate the t-t'-U model in the limit of 
t' — > in two and three dimensions and the extended 
Hubbard model in two dimensions. 

In the first part of the paper, we develop the analyti- 
cal background needed for the calculations. We end up 
with a closed set of self-consistent equations for the anti- 
ferromagnetic thermal eriuilibrium state, which respects 
the Pauli principle and the particle-hole symmetry. The 
numerical evaluation of these self-consistent equations is 
presented in the second part of the paper. The com- 
plete phase diagrams considering the paramagnetic and 
the antiferromagnetic phase for both models and in dif- 
ferent dimensionality are calculated, as well as the den- 
sity of states and the quasi-particlc dispersions. The 
main results of the paper can be considered: the pres- 
ence of three kinds of transitions (Mott-Hubbard, Mott- 
Heisenberg and Heisenberg) at half-filling in the plane 
T-U ; the existence of two components in the antiferro- 
magnetic gap (one due to the antiferromagnetic correla- 
tions and another coming from the Mott-Hubbard mech- 
anism); a finite critical value of the Coulomb interaction 
for the Mott-Hubbard and Mott- Heisenberg transitions; 
the strong decay of the Neel temperature with doping; 
a metal-insulator transition away from half-filling inside 
the antiferromagnetic phase. 



II. THE HUBBARD MODEL WITHIN THE TWO 
POLE APPROXIMATION 

The Hamiltonian for the single-band Hubbard model^ 
with chemical potential ii is defined by 

H =Y1 ~ ^^'j) (^) 0') 



+ t/^nt ii)ni (i) 



(2.1) 



with a S {Ti i} and where the sum over the site indices 
i,j runs over the whole chemical lattice {Ri}, which we 
consider to be of simple cubic type with lattice constant 
a. U is the on-site Coulomb interaction. The kinetic 
part of the Hamiltonian is given by the nearest neigh- 
bor hopping. In d dimensions the hopping matrix takes 
the form = -2dta^J = -2 dt;^ e*'' (^'-^j)a (k), 

with a (k) — ^ X^iLi cos (fc/a). 

In the framework of the COM we take as basic fields 
the Hubbard operators: ^o- (*) = Ca (i) [1 — n^- (i)] and 
Vcr (i) = Co- (i) n-a (i), where ria {i) = cj. {i) («), describ- 
ing the basic excitations n{i) = ^ n{i) = 1 and 
n{i) = 1 n{i) = 2 on the lattice site i, which are re- 
sponsible for the leading contributions in the electronic 
density of states. Let use define the following four- 
component composite field at lattice site i: 

{i,t) 

where the variable t denotes the time translate of the op- 
erator under the Heisenberg dynamics generated by H. 
The Heisenberg equation of motion for this basic field \E' 
may be split into a part that is linear in the basic field 
and the remaining non-linear part 

. d 



(2.2) 



at 



* (i, t) = j{i,t) = J2^ (i^ j) * (i' t) + *) 



(2.3) 

by projecting it on the basic field. Then, the matrix 
e{i,j) is fixed by 



e{i,j)=J2{{Jii't),^^ {l,t)}) 

x({*(M),*^ {j,t)})- 



(2.4) 



where (•) denotes the thermal average on the grand- 
canonical ensemble. 

To simplify the notational effort we introduce the so- 
called normalization matrix I and the so-called m-matrix 

= {{^ {i,t) ,^Hj,t)}) (2.5) 

m{i,j) =({J{i,t),^Hj,t)}) (2-6) 

and get the shorthand for the energy matrix e = ml~^ 
from Eq. (2.4). 

The pole approximation consists in neglecting the non- 
linear part dJ{i,t) in the equation of motion. This is 
equivalent to the suppression of the incoherent part in 
the retarded Green's function matrix 

s{i,j,t) =e{t){{^{i,t),¥{j,o)}) 



I r a 



) 



dioe 



[ d'^pe'P-^'5(k,p,a;) 



(2.7) 
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which then satisfies the linearized equation of motion 

i^^S{i,j,t)=i6{t)I{i,j) + J2^i^^^)S{l,J,t) . (2.8) 

d is the dimension of the system, a the lattice spacing 
and Qb the Brillouin zone. The resulting set of alge- 
braic equations for the Fourier transform of the retarded 
Green's fimction may then be solved. The knowledge of 
the retarded Green' function will allow us to calculate 
the correlation functions ('^ {i,t)'^^ {j,0)) by means of 
the spectral theorem. As already mentioned, the above 
scheme of truncating the equation of motion by a pro- 
jection technique is rather common in the treatment of 
systems with strong correlations. When we use as basic 
fields, as we do in this work, the operators coming from 
the hierarchy of the equations of motion,^^' the COM is 
similar to the SDA}^''^'^-'^^ However, the COM gives the 
possibility of choosing the basic fields according to the 
physics of the system under investigation. This freedom 
allows a better control on the dynamical information still 
present in the generalized mean-field approximation. Let 
us emphasize that even in the pole approximation the 
electronic self energy is not trivial as in the standard 
mean field approach, but a pole expansion of the exact 
electronic self energy is used. Another fundamental dif- 
ference between the COM and the approaches mentioned 
above consists in the treatment of the higher-order cor- 
relation functions occurring in the energy matrix (2.4), 
which are not directly related to elements of the Green's 
function. As shown in a recent publication,^*^ there is no 
freedom in choosing the equations to compute those pa- 
rameters as they have to be used to fix the representation 
according to the following relation 

lim 5(i,i,t) = (*(i)*t(i)) , (2.9) 

where the l.h.s. comes from Eq. (2.7) and the r.h.s. de- 
rives from the basic algebraic properties of the electronic 
field algebra - namely the Pauli principle. Any other 
choice of the self-consistent equations that fix those pa- 
rameters will lead to a representation for the Green's 
function where the main symmetries of the system are 
violated.^* In the case of the Hubbard model, Eq. (2.9) 
leads to the following self- consistent equations 

{U^)vt{^)) = (2.10a) 

(eT«4«) = (^iW4W)- (2.10b) 

The use of these equations, in the paramagnetic phase, 
led to a good agreement of the COM with the numeri- 
cal results for the local, integrated and thermodynamic 
quantities. ^^'^^'^^ 

In contrast to the paramagnetic solution - character- 
ized by a complete translational and spin rotational sym- 
metry -, where the number of parameters in the energy 
matrix emerging from higher-order correlation functions 



equals the number of constrains given by the Pauli prin- 
ciple (2.10), an antiferromagnetic solution with two com- 
posite fields, which is characterized by a broken transla- 
tional and spin rotational symmetry (they are, however, 
not broken in a completely independent way), requires 
additional parameters to satisfy all constrains emerging 
from the algebraic properties of the basic field opera- 
tors ^ and T]. It turns out that for the extensions of the 
simple Hubbard model (2.1) by a next nearest neighbor 
hopping term, the t-t'-U model, or by a nearest neigh- 
bor Coulomb repulsion, the so-called extended Hubbard 
model, the number of those parameters equals the num- 
ber of 'algebraic' constrains in the case of an antiferro- 
magnetic solution. According to this, in the following we 
will investigate the antiferromagnetic solution of the 2D 
and 3D t-t'-U model in the limit of t' going to zero and 
of the 2D extended Hubbard model. 

For the antiferromagnetic solution the two poles cor- 
responding to the two composite fields ^ and rj, will split 
up into four poles due to the broken translational and 
spin rotational symmetry. 

III. THE ANTIFERROMAGNETIC SOLUTION 
FOR THE HUBBARD MODEL 

A. The Antiferromagnetic Solution on a Bipartite 
Lattice 

As stated in the above we will study the antiferromag- 
netic solution of the Hubbard model by considering the 
t' limit of the t-t'-U model. The Hamiltonian of the 
latter reads as 

if«'^=^(t,^.+t;^.-/,5,^.)4 i^)c^ (j) 

^ (3.1) 
+ U2_^n^ {i)n^ {i) . 

i 

where the t'^j matrix describes the next-nearest neigh- 
bor hopping. The t'-hopping is an intra-sublattice hop- 
ping for each of the two sublatticcs A and B of a bi- 
partite square lattice, whereas the t-hopping is an inter- 
sublattice hopping. 

In this paper, wc consider an antiferromagnetic solu- 
tion of the Hamiltonian (3.1) cliaractcrized by an oppo- 
site sign in the spin density on nearest neighboring lat- 
tice sites. Thus, the antiferromagnetic ordering induces 
a magnetic lattice {Ri} with lattice constant a = a\/2, 
overlaying the chemical Bravais lattice {R-i} with lattice 
constant a. As shown in Fig. 1, the magnetic lattice is 
obtained as a square lattice with basis by collecting to- 
gether two neighboring sites of the chemical lattice. 

During the calculation of the Green's functions it will 
be convenient to switch between the two equivalent rep- 
resentations of the system constructed on the chemical 
and on the magnetic lattice. We denote vectors belong- 
ing to the chemical lattice by their bare symbols, whereas 
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the corresponding vectors on the magnetic lattice are de- 
noted by an additional '~'-symbol. 




O 



O 



O 



O 



FIG. 1. The bipartite lattice structure induced by the an- 
tiferromagnetic ordering on a 2 dimensional square lattice. 

We impose the following global boundary conditions 
on the antiferromagnetic thermal equilibrium state: 

Assumption III.l 

The antiferromagnetic tliermal equilibrium state has to 
satisfy the following global boundary conditions: 

1. The particle and spin densities are globally con- 
served. 

2. The mean value of the local particle density per 
spin satisfies the relation 

= i[n -(-!)" m cos (Q-R,)] (3.2) 

where m is the staggered magnetization, which has 
to be calculated self-consistently. The particle den- 
sity n is imposed as an external parameter. Further 
we adopt the convention (7 = 1,2 for u =t, i and 



Q 



'77 TT_) 



3. All expectation values are invariant under the con- 
temporary exchange of the spin direction and the 

sublattice index. 

4. The normalization matrix I and the m matrix from 
Eqs. (2.5) and (2.6) have to be real. 

We remark that the conditions 1.-3. are natural bound- 
ary conditions for an antiferromagnetic thermal equilib- 
rium state characterized by a staggered magnetization, 



whereas the condition 4. is a direct consequence of the 
fact that wc are looking for a state without quasi-particle 
damping, as required by the approximation scheme de- 
scribed in Sec. II. 

To reduce the computational effort in solving the in- 
tegral self-consistent equations we use a spherical ap- 
proximation for the i'-hopping, characterized by an ad- 
ditional hopping to next-next-nearest neighbors with 
half the weight of the next-nearest neighbor hopping, 
in the calculations of the integrals in the momentum 
space. This permits a definition of the t'-hopping ma- 
trix by Q:(k) only and therefore reduces the number 
of evaluations to 0{N) - N the number of steps in 
the integrals - independently from the dimensionality 
of the system, while it would have been 0{N'^) for a 
d-dimensional system with the usual next-nearest neigh- 
bor t'-hopping matrix. In 2D the t'-hopping matrix is 
given by t^^ = -12t%- = -12t'^ e''"^'^'"'^^'/^ (k) 
with /3(k) = i(4[a(k)]^ - I). In 3D we find t'^j = 
-SOt'pij = -30t'^X;ke''"^"-*"^^'^/3(k) with /3(k) = 
|(3 [a (k)]^ — i). Let us remark that in the present defi- 
nitions the hopping per lattice site is always normalized 
to 1. 

It is worth noting that, as we will take the t' ^ limit 
hereafter, the use of the spherical approximation can not 
affect at all the results we will obtain. 

In the following we will restrict our analysis to the 2D 
case. The 31? case, which might be treated in complete 
analogy by renormalizing the hopping constants and by 
changing the projections and will be considered 
in Sec. IV C. 



B. The Normalization- and the Energy- Matrices 

In order to calculate the energy matrix in the two- 
pole approximation under the antiferromagnetic bound- 
ary conditions given in Ass. III.l, we proceed along the 
guidelines given in Sec. II. To this end it is convenient to 
use the corresponding quantities defined on the chemical 
lattice and to switch to the magnetic lattice only for the 
calculation of the Green's functions. 

Firstly, we calculate the full equation of motion for the 
composite field spinor (2.2) according to Eq. (2.3) for the 
t-t'-U Hamiltonian 



/-/i^t {i, t) - Ate? {i, t) - Atn? {i, t) - 12tc^ {i, t) - m'n^ {i, t) \ 



z-^ii,t)=[^ii,t),H'*'^ 



-Urj-i {i,t)+4tTr^ {i,t) ■ 
Uc"^ {i,t)- 4t7r° {i,t) 



m'lr': {i,t) 

12tc^it,t)^ 
Ut'irf {i,t) 



m'-rrl {i,t) 



(3.3) 



4 



where we used the notation ttJ (i) = —n-a (i) (i) + 
cL{i)c^{i)cl{i) + (i) 4^ (?) (i) with cl{i) = 
'^jJijC^ij), where j = a,/3 stands for the projections 
denned in Sec. Ill A. 

According to Ass. III. 1.1, the expectation values for 
spin exchange operators vanish and we find the general 
block structure: 



(3.4) 



M = Mt©M4 , S = S^®Si. 

The explicit form of the normalization matrix is obtained 
by direct evaluation of Eq. (2.5) 

7(i,i) = + cos(Q • R,) (3.5) 

with 7(") = 7(") 7(") and 7('") = 7('") -7(™) where 
the blocks are given by 



/(") = 



j(m) _ 



2 ^ 




For the m-matrix we get from Def. (2.6) 

m(i,j) =6^JMl + ay-Ma + /Jy-Mg 

+ cos (Q • Ri) {6ijM2 + aijMi + PijMe) . 

(3.7) 

The matrices Mi , Mg also have the above mentioned 
block structure 



Ml 
M3 

M5 



Ml 


\ 


, M2= 1 


f M2 








Ml ] 




I 


-M2 


M3 


\ 


, M4= 1 


f M4 








M3 j 




i 


-M4 


Ms 


\ 




( Me 








Ms ] 




i 


-Me 



(3.8) 



(3-6) with 



Ml = 



M2 = 



Ms = 



M4 



A^) 



4t(A5' + A!j») +6t'(A^ + A^) 



-/xf - 4i(A'f - A' 



6t'{Af 



A?) 



4i(Ar 



At) 



6t{A{ 



A?) 



(i7 - /x) in - 4t(AJ + A^) - 6f' {Af + A^) 
4i(Ay-A^)+6t'(Af-A^) \ 



4t{A» - A») + 6t'{Af - Af ) -f (-M + C/) - 4t(A« - A«) - 6i'(Af - Af ) 



-4i(l-n + p) -4t(in-p) 
-4f(in-p) -4tp 



Ms = 



-12t'(l - n + i(pf +P?)) -12i'(in - i(pf + p^)) 
-m'(in-i(pf+pf)) -6t'(pf+p?) 



-2tm 
2im 



Me 



-12i'(m+i(pf-pf)) 6i'(m+pf-pf) 
6i'(m + pf-pf) -6i'(pf-pf) 



(3.9) 



To abbreviate the notation of the occurring expecta- 
tion values we defined the following parameters 

K = l{{&ii)cfi^))-{^^i^)vti^))) , for * £ A 
P=\{r^, (0 K (^)> - (^T (i) ci {i) [c\ {i) c\ (x)] ") 
A^ = (C. {i) cf (i)) - (4 {i) vl (i)) , for i e A 
= l((ni (0-(-)'^*n2(*)) 
x[ni (*) + (-)^*n2 (z)]'^') + (n. (O) 
-(c^(i)c<,(i)[4(j)4(j)]'^) forieA. (3.10) 

Furthermore, we used the notation (i) = {i) a/^c {i) 



for the spin- and charge- density operators with the Pauli 
spin matrices cr^ S {l,ax,(Jy,crz} and the electronic field 
spinors (i) ~ (c|(j) , c|(i)). 

As the antiferromagnetic ordering breaks the transla- 
tional invariance, the parameters A^, p, A^ and in 
principle do depend on the lattice site i. However, the 
antiferromagnetic state enjoys a translational invariance 
within each one of the two sublattic:es A and B. For the 
definitions of the parameters (3.10) we arbitrarily choose 
the values on the sublattice A. Their values on the sub- 
lattice B arc then given by exchanging the spin indices 
according to Ass. III. 1.3. 

The operators, from which the expectation values for 
the parameters A", p, A^ and p^ arc taken, arc not hcr- 
mitian. However, the corresponding parameters have to 
be real according to Ass. III. 1.4. This finally results in 
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the parameter p being independent of both the spin and 
the sublattice. Furthermore, the normaUzation matrix 
and the m-matrix result symmetric. 

To calculate the energy matrix we note that Eq. (2.4) 
gives m = el. In the Fourier space we have 

C^' P) = TT^ / £ (k + q, p - q) 7 (q) 

[2n) JiiB (3-11) 

^ e (k, p) + e (k + Q, p - Q) /(™) 

where the Fourier transform of the normalization and the 
m-matrix are given by 

'27r\2 



m 



/ (k, p) = j (5 (p) /(") + <5 (p - Q) /("')^ 

(k, p) = (^)' ( 5 (p) (Ml + a (k) M3 + /3 (k) Ms) 

+ ^(p-Q)(M2 + a(k)M4 + /3(k)M6)) . 

(3.12) 

We recall that the Fourier transform on the chemical lat- 
tice is defined as in Eq. (2.7) to benefit from the period- 
icity of the thermal equilibrium states. 

Using the 2Q periodicity of the thermal equilibrium 
states in Eq. (3.11) we get the energy matrix in Fourier 
space 

e (k, p) = m (k, p) C - m (k + Q, p - Q) £) (3.13) 
with 

-In -1 



Q = (/M)"^ [/(») (/M)"^ _ /M (/(")) 



In -1 



(3.14) 

Using the explicit expressions (3.12) for the m- and the 
normalization matrix we can write the energy matrix as 



e (k, p) = (1^) ' ((5 (p) + a (k) e^^^ + /3 (k) e 
+ (5 (p - Q) (e(3) + Q, ^(4) + ^ (k) e 



(5) 
(6) 

(3.15) 



with 



= MiC - M2D , e(2) = M3C + M4D 
e(3) = M2C - MiD , e(^) = M4C + M3D 
e(5) = M5C - MgD , £(^) = MgC - M5D . 

(3.16) 



C. The Green's Functions 



The solutions of the linearized equation of motion (2.8) 

with the expression (3.15) for the energy matrix can be 
interpreted as translational invariant Green's functions 
S^^(k,Lo), 5^^(k,a;), S^^i^oj) and 5^-«(k,a;) on the 
magnetic lattice as is illustrated in App. A. They have 
the general structure 



S^'^Klu) = ^^^+coA^^(k) + B^^(k) 



(a;C^^(k) + D^^(k)) 



. (3.17) 



The explicit form of the coefficients A^^ , B^^ , C^^ 
and D^^ are given in App. II. 

Due to the assumptions III.l on the antifcrromagnctic 
thermal equilibrium state, the Green's functions (3.17) 
have the block structure S^^ = Sf^ Sf^ shown in 

Eq. (3.4), where the poles {^^f (k) |i e {1,2,3,4}} of 
the spin dependent parts S^^ of the Green's functions 
are given by the roots of the fourth order equations 



det 



(a;2 + wA^^(k) + B^^(k)) = . (3.18) 



The set of poles for the Green's fimctions are all 

equal, i.e. the quasi-particle energies do depend neither 
on the spin nor on the sublattice. This reflects the prop- 
erty of the antifcrromagnctic state with staggered magne- 
tization, where the majority spin states and the minority 
spin states energetically occupy exactly the same regions 
and differ only in their corresponding spectral weights 
(cf. also Sec. IV). Therefore, we simply write for the 
poles Ei(k) = (k). 

We can write the retarded Green's functions (3.17) as 



5^^(k,a;)=limy: 

n — >u ^ * 



a^r(k) (3.19) 



^^'^^^<jj-Ei{\<i)+ir] 



with X, y e {A, B}, a € {t, i} and the spin- and 
sublattice-dependent spectral weights a^J given by 



<f(k) = ^ — 

n(i^i(k)-£;,(k)) 



{E,{^.)fc^^{^.) + (^,(k))^(D^^(k) + det(^^^(k)) (A^^(k))-^C,^^(k)) 



^ + i?.(k)(det(A^^(k)) {A^^{k)y'Dr{k) 



+ det(B--(k))(B--(k))-^C--(k) + det(i?--(k))(B--(k))-^2)--(k) 



(3.20) 
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D. The Self Consistency Equations 

Given the temperature T and the particle density n 
as external thermodynamic parameters as well as the 
Coulomb interaction U and the hopping constants t and 
t' , which are the model dependent parameters, we are 
now able to give a closed set of self-consistent condi- 
tions for the internal parameters characterizing an an- 
tiferromagnetic thermal equilibrium state. These inter- 
nal parameters are p, p^, p^, A", A", and A^ from 
Eq. (3.10), the chemical potential fj. and the magnetiza- 
tion m from Eq. (3.2). 

For the calculations of these parameters we need the 
knowledge of the correlation functions, which are con- 
nected to the retarded Green's functions of the funda- 
mental spinor by means of the spectral theorem. In view 
of the special form (3.19) of the retarded Green's func- 
tions the spectral theorem at equal time may be written 
as 



C^^(Ri,R,-) = (*^(RO*^+(R,) 



J2 / d'ke^<^'-^' 



X a, 



XY 



(k) 



1 + tanhl 



(3.21) 



We denote the on-site, the nearest neighbor and the next- 
(next)-nearest neighbor correlation functions at equal 

time by C^^(Ri), C^^"(Ri) and C^^^(RO. 

For the parameters /U and m we then find the self- 
consistent equations 



2-n==CiY(R.)+C^n''(R.) 
2m = Qy(Ri)-C4-4(Ri) 



(jAA 



CP"" 



22 (f^i) 
22 (R-i) 



C^^(Rj)] . 



22 (R-i) 
44 



(3.22a) 



The parameters A", AJ^, and A^ are directly related 
to matrix elements of the correlation functions by 



A? = 



A? = 



A^ 
Af 



[C^\^"(Ri) + Ci^2^"(R0 

JC3i^«(R,) + C3^/"(Ri) 

— O34 (±tij - O44 
[C,r^(R,)-C2l^^(R,)] 

[C^^'\^i) - C^'^'^^(Ri)] . 



(3.22b) 



The parameters p, p^ and p^^ cannot be calculated ex- 
plicitly by the single-particle Green's function (3.19), be- 
cause they derive from higher-order correlation functions. 



According to what stated in Sec. II, we will use the fol- 
lowing equations to fix the representation of the Green's 
function^^'^^ 



(7i^2^(R0 = 



(3.22c) 



Cf2^(Ri) = 

We observe that all the self-consistent equations are cou- 
pled and have to be solved as one set by means of a global 
convergency scheme. 

Finally we remark that the correlation functions 
C^^(Ri), C^^^(Ri) and C^^^(Ri) actually do not de- 
pend on the lattice site Rj of the magnetic lattice because 
of the translational invariance enjoined by the Green's 
functions (3.19). 



IV. NUMERICAL EVALUATION OF THE 
ANTIFERROMAGNETIC PHASE FOR THE 
HUBBARD MODEL 

In the following the thermodynamics of the antifer- 
romagnetic thermal equilibrium states for the Hubbard 
model in two and three dimensions are discussed. We 
report the phase diagrams resulting from solutions of 
Eqs. (3.22) as well as some of the microscopic properties 
of the corresponding solutions, which explain the phase 
behavior. In this context, our main concern is in the in- 
terplay of the antiferromagnetic (Mott-Heisenberg) gap 
and the Mott-Hubbard gap leading to a metal-insulator 
transition in the antiferromagnetic phase. Furthermore, 
we investigate the distribution of spectral weight between 
the majority and minority spin states that explains most 
of the phase properties found for the antiferromagnetic 
solution. 

As a first step, we establish an averaging procedure be- 
tween solutions for positive and negative values of t' that 
is capable to rule out a non physical artifact induced by 
the spherical approximation and to establish the corre- 
spondence to the simple Hubbard model. 



A. The Averaging Procedure in t' 

Due to the incompatibility of the nearest and the next- 
nearest neighbor hopping the t-t'-U model does not en- 
joy the particle-hole symmetry at half filling {n = 1). In 
addition, the spherical approximation for the t'-hopping 
overemphasizes the intra-sublattice hopping to the next- 
and next-next-nearest neighbors. This leads to an in- 
stability of the antiferromagnetic solution for the t-t'-U 
model at half filling: the magnetization goes to zero for 
positive values of the f'-hopping while for negative values 
of t' diverges to infinity (cf. Fig. 2(a)). Positive values 
of the intra-sublattice hopping do suppress antiferromag- 
netism to a certain extent, whereas negative values are 
in favor of it by the additional phase factor of tt. 
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For the simple Hubbard model the particle-hole symme- 
try reflects in the algebraic property 

M(n = 1) = I . (4.1a) 

For an antiferromagnetic solution in the framework of the 
two-pole approximation we have the additional condition 

A«(n = 1) = -A'^in = 1) , (4.1b) 

which is a generalization of the condition A"(n = 1) = 
= A"(n = 1), imposed by the particle- hole symme- 
try on the paramagnetic solution of the simple Hubbard 
model, because of the inequivalence of the majority and 
the minority spin subsystems in the antiferromagnetic 
state. The condition (4.1a) is satisfied in the limit t' go- 
ing to zero independently of the direction (see Fig. 2(b)). 
This is a direct consequence of the fact that the proper 
representation for the Green's functions has been taken 



by using the complete set^°' of constrains coming form 
the Pauli principle (3.22c). 

Taking the average between the solutions for t' and —t' 
with \t'\ approaching zero we get an antiferromagnetic 
state - which formally has zero t'-hopping - that satis- 
fies Eqs. (4.1) and has no divergence in the magnetization 
and the other parameters at half filling (see Figs. 2(a), 
2(c) and 2(d)) and can be considered as representative 
for the simple Hubbard model. Averaging the solutions 
for ±t' thus combines the large benefit in computational 
time provided by the spherical approximation with an 
antiferromagnetic solution that satisfies the complete set 
of symmetry constrains deriving from the Pauli principle 
and enjoys the particle-hole symmetry. The numerically 
accessible limit for the t'-hopping, which has been used 
in combination with the above averaging procedure (see 
Fig. 2), is t' = ±10~^i. Hereafter, we will present re- 
sults exclusively from the averaged solution and we will 
consider them as those of the simple Hubbard model. 
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FIG. 2. The averaging algorithm between positive and negative values oft' in the 2D t-t'-U model for U = lOt and kT = 0.5t 



8 



B. The Antiferromagnetic State of the Hubbard 
Model 

1. The Phase Diagram 

The n-T and U-n phase diagrams for the antiferro- 
magnetic state of the 2D Hubbard model and the cor- 
responding paramagnetic state are shown in Figs. 3 and 
4. The antiferromagnetic state has a free energy lower 
than the one of the paramagnetic state over the whole 
phase region, leading to a phase transition of second or- 
der between the antiferromagnetic and the paramagnetic 
phase at the lines of vanishing magnetization. The study 
of the phase diagram near half filling has to be completed 
by an investigation of the ferromagnetic phase^° and the 
charge ordered phase, which both can be studied in the 
framework of the approximation scheme described above. 
The antiferromagnetic phase could be energetically ruled 
out by one of these other phases in certain regions of the 

phase diagrams, then leading to phase transitions of first 
order.26. 31-33 

In this section we give a brief overview of the proper- 
ties of the antiferromagnetic phase, which then will be 
related to the inner structure of the antiferromagnetic 
state - namely its density of states - in the subsequent 
sections. 

The most striking features of the antiferromagnetic 
phase are the finite critical Coulomb interaction Uc as 
a lower bound to the antiferromagnetic state at half fill- 
ing, the restriction of the antiferromagnetic state to a 
very narrow region in n around half filling and the metal- 
insulator transition (MIT) within the antiferromagnetic 
phase. The vanishing of the staggered magnetization at 
half filling for U < Uc - we find Uc within 5t-10t (see 
Fig. 4) - is supposed to be an effect of strong electron cor- 
relations. It cannot be observed within simple mean- field 
treatments of the Hubbard model, ^^"^^ where the antifer- 
romagnetic phase is stable down to t/ = at n = 1. Also 
in some more sophisticated mean-field approximations, 
as the SDA, a stable antiferromagnetic state is found at 
n = I down to very small values of U and its stability 
down toU = cannot be excluded. The same holds for 
the antiferromagnetic state of the simple Hubbard model 
treated by the DMFT.^'^ 

As we already mentioned in Sec. I, the SDA is rather 
closely related to the COM, but to calculate the anti- 
ferromagnetic state further approximations on the corre- 
lation functions were needed. ^^'^^ The main difference 
lies in the treatment of the internal parameters emerg- 
ing from higher-order correlation functions. While the 
COM uses them to fix the representation of the Green's 
functions (3.22c), they are calculated by the equations of 
motion in the SDA. 

In a previous study^^ we investigated the critical 
Coulomb interaction as a function of the t'-hopping, 
showing that a finite value of Uc in the above mentioned 
range has also to be expected in the exact limit t' — > 0. 




0.94 0.96 0.96 0.97 0.98 0.99 1 0.94 0.95 0.96 0.97 0.98 0.99 1 

n n 
(a) U = lot (b) U = 20t 

FIG. 3. The n-T Phase diagram for the 2D Hubbard 
model. 




(a) kT = O.Olt (b) kT = 0.5 

FIG. 4. The Phase diagram for the 2D Hubbard 

model. 



The antiferromagnetic phase found in the high- Tc cop- 
per oxide compounds in general shows a great stability 
at half filling, whereas already a few percent of electron 
or hole doping leads to a strong reduction of the Neel 
temperature and eventually to the vanishing of the anti- 
ferromagnetic phase. On the other hand it is well known 
that in general the mean-field treatment of the Hubbard 
model strongly overemphasizes the stability of the an- 
tiferromagnetic phase in doping. In the COM, where 
the strong correlation effects are restored to some extent 
by the Pauli principle symmetry constraint, the stabil- 
ity of the antiferromagnetic phase is actually reduced 
to a narrow region of a few percent of doping around 
half filling with a strong reduction of the Neel tempera- 
ture. For the Hubbard model such a behavior has been 
confirmed by numerical results in quantum Monte Carlo 
studies.'^®' ^® In contrast a larger region for the antiferro- 
magnetic phase has been obtained by the SDA^^ and the 
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d = oo approximation.^'"' 

At zero temperature a transition from an antiferro- 
magnetic insulator at half filling and an antiferromag- 
netic metal for n < 1 is observed. At higher temper- 
atures, however, we find an extended region in doping 
around n = 1, where the Fermi level is situated inside 
the Mott-Heisenberg gap, which itself is large with re- 
spect to the thermal energy kT. We thus have a state 
with poor conductivity of semiconductor type, which we 
call the 'antiferromagnetic insulator'. The metallic be- 
havior is recovered when the Fermi level joins the second 
antiferromagnetic band with decreasing values of n. 

The n-T phase diagrams show a strong qualitative dif- 
ference in the low temperature region according to the 
presence (Fig. 3(b)) or the absence (Figs. 3(a)) of a Mott- 
Hubbard gap. For values of U that do not admit a Mott- 
Hubbard gap the stability of the antiferromagnetic phase 
is enhanced by increasing temperature (a phenomenon 
called 'heat magnetization'), whereas the contrary is true 
for values of U where the Mott-Hubbard gap is already 
opened. 

Furthermore, the presence of the Mott-Hubbard gap 
and the fact that it closes within the antiferromagnetic 
phase, when the particle density is reduced, is also re- 
sponsible for the intermediate reduction of the stability 
of the antiferromagnetic state by increasing the Coulomb 
interaction (Fig. 4). This reduction is reinforced by in- 
creasing temperature and leads to a 'nose-like' shape of 
the U n phase diagram. Both phenomena can be ex- 
plained by the evolution of the spectral weights for the 
majority and minority spin subsystems as functions of 
the external parameters, as will be explained in detail in 
the following section. 

Finally, we remark that the phase diagrams shown in 
Figs. 3 and 4 are completely symmetric with respect to 
n — 1] this is due to the fact that the thermal equilibrium 
state respects the particle-hole symmetry. To incorporate 
the experimentally observed asymmetry in particle- and 
hole-doping a projection of the two-band Hubbard model 
on an effective single-band one has been proposed. 



2. The Band Properties 

The spectral properties of the antiferromagnetic state 
arc deduced from the electronic single-particle Green's 
function on each sublattice S^^^ {k, uj) with X e {A, B}. 
In the following we will restrict our analysis to the sub- 
lattice A, while the quantities on the sublattice B are 
obtained by simply exchanging the majority and minor- 
ity spin subsystems. We thus have for the retarded elec- 
tronic Green's function on the sublattice A 

SifiKu;) =5fjf (k,a;) + Sf^^{k,c.J 



The spin dependent electronic density of states within 
the antiferromagnetic state is then given by 

"'<-'^=(l^)'4''*(-^'"'*"<*'^') 

•'^'s i=l 

X f<^,t(k,a;)+a1^/ t(k,t^) 

+ a^^,t(k,a;)+a^^ t(k,a')) ■ 

(4.3) 

We recall that all energies are referred to the chemical 
potential. 

The two-pole approximation leads to a splitting of the 
single electronic band into two Hubbard subbands that 
correspond to the elementary excitations described by 
the composite operators ^ and 77. Those two subbands are 
separated by the on-site Coulomb interaction [/, which 
may lead to a gap at same critical value. The bi- 
partite lattice approach leads to a doubling of the two 
subbands by reflection around the band center at — A* 
(see Fig. 6(c)). In the paramagnetic case, when the mag- 
netic lattice is introduced, this doubling of the Hubbard 
subbands by the reduction of the Brillouin zone is com- 
pletely artificial: the reflected bands occupy exactly the 
same energy interval with the same spectral weights as 
the original Hubbard subbands. 

The antiferromagnetic state, however, is characterized 
by the opening of Mott-Heisenberg gaps at the crossing 
points of the four subbands induced by the switching to 
the magnetic lattice. Thus we find three Mott-Heisenberg 
gaps, one in the lower Hubbard band, A^, one in the up- 
per Hubbard band, and a central Mott-Heisenberg 
gap, Aj_^, in the region where the two Hubbard sub- 
bands overlap in the paramagnetic phase (see Fig. 5). 
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FIG. 5. The density of states on the sublattice A, illus- 
trating the Mott-Heisenberg gaps in the lower and the up- 
per Hubbard band as well as in their overlapping region in 
the paramagnetic phase. The values of the parameters are: 
U = lot, n = 1.0 and kT = O.Olt. 
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If the two Hubbard subbands arc already separated 
by a Mott-Hubbard gap (in the paramagnetic phase it 
has been found^" that the Mott-Hubbard gap opens at 
U ~ 13. 2i), the eentral Mott-Heisenberg gap adds to this 
gap (see Fig. 6(a)). The gaps and A^, are not sym- 
metric around the Mott separation ±[//2, because the 
upper and the lower band edges are shifted by different 
amounts. 

The antiferromagnetic ordering thus leads to a split- 
ting of the two Hubbard subbands into four antiferro- 
magnetic bands that are occupied by both, the majority 
and the minority spin subsystems. As it can be seen from 
Fig. 5, the two spin subsystems energetically occupy the 
same regions - which is the reason why the staggered 
magnetization does never reach saturation m = n -, but 
with rather different spectral weights: the spectral weight 
of the majority spins is strongly enhanced at the upper 
band edges, whereas the minority spins have an enhanced 
spectral weight at the lower band edges, leading to the 
staggered magnetization. At half filling the density of 



states shows a complete symmetry between the minority 
and majority spins with respect to reflection of the en- 
ergy around the Fermi level. This is due to the fact that 
our solution respects the particle-hole symmetry. 

a. The Mott-Heisenberg and the Mott-Hubbard gap 
The interplay between the Mott-Hubbard gap and the 
Mott-Heisenberg gap at half filling is illustrated in Fig. 6. 
For U = lot the two Hubbard subbands overlap at the 
Nccl temperature Tn and by decreasing T we find the 
opening of the Mott-Heisenberg gaps in the two Hubbard 
subbands as well as in the central region (Fig. 6(b)). The 
corresponding evolution of the electronic band structure 
in the Brillouin zone of the chemical lattice is shown in 
Fig. 6(c). We notice the doubling of the Hubbard sub- 
bands in the paramagnetic state, which then evolves into 
the four antiferromagnetic subbands by decreasing tem- 
perature, and the typical symmetry of the band structure 
along the diagonal of the Brillouin zone for the chemical 
lattice due to the reduced Brillouin zone of the magnetic 
lattice. 
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(a) Sublattice Density of States for 

U = l&t 



(b) Sublattice Density of States for 

U = Wt 



(c) Band Structure on the Chemical 
Lattice for U = Wt 



FIG. 6. Density of states and band structure at half filling; the thin vertical line at B = denotes the position of the Fermi 
level Ef- 



In Fig. 6(a) the case for U = l&t is shown, where 
the Mott-Hubbard gap already separates the two Hub- 
bard subbands in the paramagnetic phase just above the 
transition point. By decreasing temperature the Mott- 
Heisenberg gaps A^ and A^ open within each of the 
two Hubbard subbands. In addition the central Mott- 
Heisenberg gap adds to the Mott-Hubbard gap leading 
to the central gap Aj_^. 

The temperature dependence of these gaps is illus- 
trated in Figs. 7(a) and 7(b) for various values of U . 
Above a certain critical value of U the central gap re- 
mains open at the phase transition. At low temperatures 



the Mott-Heisenberg gaps within the two Hubbard sub- 
bands, A J and A^, are increasing as a function of U up 
to a maximum value U ^ lit and then arc decreasing like 
1/U (cf. Fig. 7(c)). This reflects the l/J7-dependence of 
the antiferromagnetic exchange integral in the Heisen- 
bcrg model, to which the half filled Hubbard model can 
be mapped. ^^'^^ For the central gap the additional Mott 
separation prevents this behavior. Also the Neel temper- 
ature does not show an 1/U like behavior (cf. Fig. 8), 
which indicates that the antiferromagnetic exchange in- 
tegral is not always properly taken into account. 
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(c) Gaps in the Lower and the Upper Hubbard Band, 
A^ and A,,, as Function of U 

FIG. 7. The Mott and the Mott-Heisenberg gap at n = 1. 



b. The Metal Insulator Transition The phenomenon 
of metal insulator transitions has been intensively studied 
in the context of strongly correlated electron systems.''"^ 
In the framework of the COM the MIT within the para- 
magnetic state at half filling is found^° in the light of the 
Hubbard picture, due to separation of the two Hubbard 
subbands at a critical value of the Coulomb interaction. 
This picture is quite different from the one found by the 
DMFT, where the MIT from a paramagnetic metal to a 
paramagnetic insulator is mainly due to the vanishing of 



a narrow coherent quasi-particle peak at the Fermi level 
as in the Gutzwiller approximation.^ However, for an an- 
tiferromagnetic state on a bipartite lattice, this interpre- 
tation can be maintained within the DMFT only when 
frustration by an additional t'-hopping is introduced.^' "^^ 
At half filling, within the framework of the COM, we 
find three kind of transitions: a Mott-Heisenberg tran- 
sition (i.e., a transition between a paramagnetic metal 
and an antiferromagnetic insulator) mainly driven by the 
temperature at low values of the Coulomb interaction; 
a Mott-Hubbard transition (i.e., a transition between a 
paramagnetic metal and a paramagnetic insulator), al- 
most insensitive to the temperature, at high values of 
the Coulomb interaction; an Heisenberg transition (i.e., 
a transition between an antiferromagnet and a paramag- 
net), within the insulating phase, driven by the tempera- 
ture at high values of the Coulomb interaction. Moreover, 
the central gap has two components depending on the 
value of the Coulomb interaction and the temperature: 
one due to the antiferromagnetic correlations (in the anti- 
ferromagnetic insulating phase) and another coming from 
the Mott-Hubbard mechanism (in the paramagnetic and 
antiferromagnetic insulating phases). In the Heisenberg 
transition the antiferromagnetic component of the cen- 
tral gap vanishes (i.e., the magnetization disappears and 
the lateral gaps close up), but the paramagnetic compo- 
nent remains finite. We have a finite critical value of the 
Coulomb interaction for the Mott-Hubbard transition, 
in contrast to what found, for instance, by the Hubbard 
I approximation and the SDA. This fact allowed us to 
study the Mott-Heisenberg transition existing at lower 
value of the Coulomb interaction and that is obviously 
absent in any picture based on the approximations men- 
tioned above (they do not have the Mott-Hubbard tran- 
sition neither!). In Fig. 8 we summarize the transitions 
occurring at half filling within a treatment of the Hub- 
bard model in the framework of the COM. 
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FIG. 8. The variety of transitions at n = 1 
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The filling controlled MIT is discontinuous at zero tem- 
perature. We have an antiferromagnetic insulator at 
n = 1 and an antiferromagnetic metal at n < 1, be- 
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cause the central gap is compensated by the disconti- 
nuity in the chemical potential at n = 1 and thus the 
Fermi level always lies inside the second antiferromag- 
netic band for n < 1. This is exactly the same result as 
found in the usual mean- field approximation'*^ and in the 
DMFT without frustration. ^^In the frustrated case with 
non-zero f'-hopping DMFT as well as Quantum Monte 
Carlo studies additionally led to a.U controlled MIT in- 
side the antiferromagnctic phase at n = 1.44,45 

Moving to finite temperature, the variation in the 
chemical potential, and thus in the overall band shift, is 
considerably moderated for higher values of T . whereas 
the central gap still remains large in comparison to kT. 
Therefore, we find a finite region around half filling, 
where the Fermi level is still situated inside A^_^ and 
we thus have an antiferromagnetic phase of semiconduc- 
tor type with very poor conductivity. When the Fermi 
level is crossing the peak in the density of states at the 
upper edge of the second antiferromagnetic band, a huge 
jump in the number of carriers is observed and we finally 
get an antiferromagnetic metal (see Fig. 9). Note that 
for the purely paramagnetic MIT such a filling controlled 
transition is not to be expected. 
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(a) kT = 0.2t (b) kT = 0.5t 

FIG. 9. The sublattice density of states a,t U = lOt for var- 
ious values of the particle density. The thin vertical line at 
-E = denotes the position of the Fermi level Ep. 

Comparing Figs. 9(a) and 9(b) we see how the mod- 
erate evolution of the chemical potential as function of 
n leads to a larger extension of the insulating antiferro- 
magnctic phase at higher temperatures. 

The almost linear dependence of the central gap Aj_^ 
on the Coulomb interaction U (cf. Fig. 7(a)) leads to a 
strong stabilization of the insulating antiferromagnetic 
phase at n = 1. Near half filling, however, this effect 
is compensated for larger values of U due to the strong 
decay of the chemical potential when moving to lower 
particle densities. This decay of the chemical potential 
then leads to a large overall shift of the band structure 
towards higher energies, bringing the Fermi level close to 
the upper edge of the second antiferromagnetic band. Af- 
ter a small region, where the insulating phase is growing 



with t/, this effects restricts the insulating antiferromag- 
netic phase to regions very close to half filling for large 
values of U. 

c. The Shape of the Phase Diagram For large 
values of the Coulomb interaction close to the critical 
value Uc, where antiferromagnetism is vanishing, the re- 
gion of filling where the antiferromagnetic phase exists is 
enlarged by increasing U. On the contrary, for interme- 
diate values of U the antiferromagnetic region in doping 
shrinks. This is can be explained by the closing of the 
Mott-Hubbard gap as a function of doping. When at low 
values of U no Mott-Hubbard gap is present, the magne- 
tization and the extension of the antiferromagnetic phase 
in n grow with growing U. After the opening of the Mott- 
Hubbard gap we find the system in a region where the 
Mott-Hubbard gap closes in the proximity of n = 1, i.e., 
already within the antiferromagnetic phase. This leads 
to a shift of the Fermi level towards the middle of the sec- 
ond antiferromagnetic band where majority and minority 
spins nearly have the same spectral weight and thus to 
a considerable reduction of the staggered magnetization 
(cf. Fig. 10). For large values of U the Mott-Hubbard gap 
remains open within the whole region of doping, where 
the antiferromagnetic phase exists, and we again find in- 
creasing stability of the antiferromagnetic phase with in- 
creasing Coulomb interaction. 
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FIG. 10. The sublattice magnetization as function of n at 
kT = 0.2t. 

d. The Heat Magnetization For small values of U we 
find an increasing of the stability in n of the antiferromag- 
netic phase by increasing temperature. In the same man- 
ner, when n is close to the phase boundary, the magneti- 
zation itself first increases with increasing T before going 
down to zero, as is shown in Fig. 11. This so-called 'heat 
magnetization' has its explanation in a strong change in 
the spectral weights of the spin subsystems by increasing 
temperature, whereas the Mott-Heisenberg gaps and the 
chemical potential are only subject to very little changes 
(see Fig. 12). At higher values of U, after the opening 
of the Mott-Hubbard gap, the stronger Coulomb interac- 
tion considerably stabilizes the antiferromagnetic state 
at low temperatures and excludes the effect of heat mag- 
netization. 
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FIG. 11. The sublattice magnetization as function of T and 
na,tU= lot. 
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C. The Antiferromagnetic State of the 3D Hubbard 
Model 



The antiferromagnetic phase of the 3D Hubbard model 
is very similar to the one observed in 2D. The shape and 
main features of the phase diagrams (Figs. 13 and 14) 
remain unchanged (i.e., the MIT within the antiferro- 
magnetic phase, the finite critical Coulomb interaction 
for the vanishing of the antiferromagnetic phase at half 
filling, the favoring of the insulating state at high temper- 
atures and low values of the Coulomb interaction and the 
heat magnetization). Furthermore, as in the 2D case, all 
these properties can be explained by analyzing the den- 
sity of states and the energy spectra. 

The higher coordination number of the 3D model re- 
duces the fluctuations and leads to a greater stability of 
the antiferromagnetic phase as a function of the external 
parameters T, n and U. Anyway, comparing the exten- 
sion in n of our antiferromagnetic phase with the one 
found within the SDA for a fcc-lattice^® the former is 
restricted to a much smaller region in doping. 



(a) Physical Phase Dia- 
gram at kT = O.Olt 



(b) Complete Phase Dia- 
gram at kT = 0.5t 



FIG. 14. The phase diagram for the 3D Hubbard 

model 



The greater stability of the antiferromagnetic phase 
shows, within the range of physically relevant values of 
the parameter U, that the metallic antiferromagnetic 
phase extends to much lower values of n (cf. Fig. 13(b) 
and 14(b)) than it might be expected from the phase di- 
agrams in Figs. 13(a) and 14(a). For higher values of the 
Coulomb interaction, up to J7 = 30t, this additional an- 
tiferromagnetic region is separated by the antiferromag- 
netic phase near half filling by a paramagnetic region. 

The form of the sublattice magnetization as a function 
of n for different values of the Coulomb interaction is 
given in Fig. 15. This result indicates that the two anti- 
ferromagnetic regions in the complete phase diagrams of 
Figs. 13 and 14 arc actually two parts of a single antifer- 
romagnetic solution, which are only connected at rather 
high values of U. We will discuss in more details the 
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nature of this 'tail' in the next section, in the context of 
the analysis of the extended Hubbard model. Here we 
only remark that the separation of the two antiferromag- 
netic phase regions by a paramagnetic one is due to the 
closing of the Mott-Hubbard gap near n = 1. This leads 
to a strong suppression of magnetization, as described in 
Sec. IVB 2c for the 2D case. The 'antiferromagnetic tail', 
which can be considered as an artifact of the employed 
approximation, is also present in the 2D case, but only 
appears at values of the Coulomb interaction U ~ lOOt. 
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FIG. 15. The sublattice magnetization as function of n at 
kT = 0.5t for different values of the Coulomb interaction. 



As we already mentioned in Sec.IVB2, this 'tail' at 
lower particle densities could be energetically ruled out 
by other phases with magnetic or spatial ordering. 

Finally, in Fig. 16, we show the evolution of the sub- 
lattice density of states by changing the Coulomb inter- 
action. We note a broader structure in comparison to 
the 2D model. As shown in Fig. 16(a), the transitions 
from metallic to sciniiconductor type antiferromagnetic 
states are related to the interplay between the global 
band shift and the evolution of the central gap. This ex- 
plains the form of the insulating antiferromagnetic phase 
in Fig. 13(a). 
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(a) n = 0.9 (b) n = 1.0 

FIG. 16. The evolution of the sublattice density of states 
under variation of the Coulomb interaction at kT = O.Oli. 
The thin vertical line at E = denotes the position of the 
Fermi level Ep. 



V. THE EXTENDED HUBBARD MODEL 

A. The Model 

In this section we study the antiferromagnetic solution 
of the extended Hubbard model in two dimensions. The 
model is described by the following Hamiltonian 

= i^ij - l^^ii) 4 (0 {j) + C/ E ii) «i (^) 

ij i 
a 

+ YVtjn„ {i)n„' (j) 

ij (5.1) 

with the nearest neighbor Coulomb interaction 

V,j = 2Vaij = 2^1 Ee^'' (^-^^)a (k) . (5.2) 

k 

We proceed along the guidelines given in Sec. Ill to ob- 
tain the antiferromagnetic state. The equation of motion 
for the fundamental spinor (2.2) is given by 



- 4te' 



^ {i, t) - Atnf {t, t) + 4T/n" [i, t) [i, t) \ 
Ur]i {i, t) + 4:tnf (i, t) + 4Vn°' {i, t) r}^ {i, t) 



—l-oill -r u //-f yi, l) t it/ii yi, l) ^ 'iv rr {i, t) r]-\ {i, t) 
-fi^l {i, t) - 4tcJ {i, t) - At-Kl {i, t) + AVn" {i, t) £,i {i, t) 
\ -mi {i, t) + UTJi {i, t) + 4t7rf {i, t) + 4Fn« {i, t) rji {i, t) ) 



(5.3) 



where n" {i) = Ylja "^y ^-ir ij)- The normalization matrix which are described by the projections (iij. The matrices 
/ takes exactly the form given in Eqs. (3.4) and (3.6). 
For the m-matrix we get the same result as given in 
Eq. (3.7) without the terms coming from the t'-hopping, 



-Ml, . . . , M4 are built up from matrices Mi, . . . , M4, as 
in Eq. (3.8), and these latter have the following form 
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Ml 



Ma 



-m(1 - ^n) - 4i(A« + A") + W{n - - Vf) 
4t(A« + A«) 

-/xf -4t(AJ-Af)+4y(-yi + yr) 
4i(A«-A«) 



4i(A? + A«) 
([/ - /x) in - 4t (A« + A«) + 4y (t;^ + v^) 



-f (-M + [/) - 4i(A« - A«) + 4y(t;^ - i;^) 



M3 = 



-4t(l - n + p) + 4yA -4t(in - p) + 4F(i^| + 1^4) 



-4ip + 4^^ 



(5.4) 



Ma 



-2tm + 4V{vi - lyi) 



where we define the additional parameters 



for I G A 



for i e A 



for z e A 



for i € A . 



(5.5) 

We remark that the parameters A and k are not spin 
dependent because of the requirement for the m-matrix 
to be real and the additional symmetry constraint in 
Ass. III. 1.3. Furthermore, the parameters A" can be 
expressed through the parameters from Eq. (5.5) as 



(5.6) 



The single-particle retarded Green's functions are cal- 
culated as in Eqs. (3.17) and (3.19), just omitting the 
parts corresponding to the t'-hopping. The parameters 
in Eq. (5.4) can then be calculated self-consistently by 
means of the correlation functions (3.21). The parame- 
ters m and /i arc calculated as in Eq. (3.22a). The param- 
eters A, i^a and k arc directly related to the single-particle 
retarded Green's functions by 



A 



''I 



1 



BAa 
12 



(5.7) 



K = C. 



ABa/ 



^22 (R-j) 

and the parameters p and Va, which derive from higher- 
order correlation fimctions arc used as in Eq. (3.22c) to 
satisfy the algebraic relations corresponding to the Pauli 
principle on the level of thermal equilibrium states 

Cil^(RO = 



C^^(Ri) 

■lAA 







(R-i) — (R; 



Again, all the self- consistent equations are coupled. 



The Antiferromagnetic State of the 2D Extended 
Hubbard Model 



The phase diagram corresponding to the antiferromag- 
netic thermal equilibrium state obtained as solution of 
the self-consistent equations (5.7), (3.22a) and (3.22c) is 
shown in Fig. 17 where we took a fixed ratio U/V = 5. 
Again, the antiferromagnetic phase always has a lower 
free energy with respect to the paramagnetic phase and 
the phase transition is always of second order. 

The antiferromagnetic thermal equilibrium state ful- 
fills the particle-hole symmetry. In addition, the main 
qualitative features of the antiferromagnetic states for 
the simple Hubbard model in two and three dimensions 
can be found in the extended Hubbard model, too. There 
is a critical value of the Coulomb interaction, for which 
the antiferromagnetic phase disappears and a metal in- 
sulator transition occurs. 

The phenomenon of heat magnetization is more pro- 
nounced than in the simple Hubbard model and the oc- 
currence of a 'tail' in the antiferromagnetic solution is 
observed down to rather low values of U, where they can 
even join for higher temperatures. This leads to para- 
magnetic inclusions within the antiferromagnetic phase. 
For higher values of U the antiferromagnetic phase is re- 
duced to a very narrow region near half filling. 

The band structure turns out to be somehow differ- 
ent from the one found for the simple Hubbard model. 
The presence of the inter-site Coulomb interaction leads 
to a large overlap of the antiferromagnetic bands, such 
that the opening of the Mott-Heisenberg gap at the cross- 
ing points cannot split the lower and the upper Hubbard 
band in most part of the antiferromagnetic phase region. 
The sublattice density of states for the extended Hubbard 
model is thus mainly characterized by a single central 
gap, which is the superposition of the Mott-Heisenberg 
and the Mott-Hubbard gap (cf. Fig. 18). 
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FIG. 17. The phase diagrams for the extended Hubbard model with V = U/5. 



The part of the antiferromagnetic phase at lower parti- 
cle density is always characterized by this overlap of the 
antiferromagnetic subbands emerging from the lower and 
the upper Hubbard band respectively. The absence of 
the Mott-Heisenberg gap within the Hubbard subbands 
reduces the number of peaks in the density of states to 
four: 




(a) Sublattice Density of 
States for U = 9M 



(0,0) (71,71) (71,0) (0,0) 

k 

(b) Band Structure on the 
Chemical Lattice for U = 
9.5t 



FIG. 18. Sublattice density of states and band structure at 

kT = OAt and with V = U/5; the thin vertical line at E = 
denotes the position of the Fermi level Ef- 



Two peaks are due to the central gap and two peaks are 
reminiscent of the van Hove singularity in each Hubbard 
band^*5 (cf. Fig. 18). 
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FIG. 19. Sublattice magnetization for different values of 
the Coulomb interaction at kT = OAt with V = U/5. 



In Fig. 19 we plotted the sublattice magnetization as a 
function of the particle density for various values of the 
Coulomb interaction U showing the separation of the an- 
tiferromagnetic phase by a paramagnetic region for high 
as well as for low values of U. 



The inter-site Coiilomb interaction V considerably re- 
duces the magnetization and the stability of the antifer- 
romagnetic state, i.e., its extension in n as can be seen 
from Fig. 20. The reduction of the magnetization also 
suppresses the separation of the two parts of the anti- 
ferromagnetic states. However, the suppression of anti- 
ferromagnctism by the inter-site Coulomb interaction is 
much smaller than the one found in the treatment of the 
simple Hubbard model. 
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FIG. 20. Sublattice magnetization as function of the par- 
ticle density for different values of the inter-site Coulomb in- 
teraction V, with f/ = 9t and kT = 0.4t 



VI. CONCLUSIONS 

Among the variety of analytical methods, developed 
in the last decades to deal with strongly correlated 
electron systems, the COM has been rather success- 
ful in describing the properties of many correlated sys- 
tems.i^'i^'19'28.40,47 a,dd another piece to the puz- 
zle constitiited by the phase diagram of the Hubbard 
model, we here investigated the antiferromagnetic phase 
characterized by a staggered magnetization. A fully self- 
consistent treatment, respecting the symmetry constrains 
emerging from the Pauli principle, has been presented for 
the Hubbard model in two and three dimensions and the 
two-dimensional extended Hubbard model. 

The antiferromagnetic phases of the three systems, 
where compared to the corresponding paramagnetic 
phases, lead to a phase transition of the second order. 
In all cases the antiferromagnetic states have lower free 
energy. Furthermore, the antiferromagnetic states of all 
the models show the following properties: 

(a) the presence of three kinds of transitions (Mott- 
Hubbard, Mott-Heisenberg and Heisenberg) at 
half-filling in the plane T-U\ 

(b) the existence of two components in the antiferro- 
magnetic gap (one due to the antiferromagnetic 
correlations and another coming from the Mott- 
Hubbard mechanism); 

(c) a finite critical value of the Coulomb interaction 
for the Mott-Hubbard and Mott-Heisenberg tran- 
sitions; 

(d) the antiferromagnetic phase is stable only in a very 

narrow region around half filling, showing a strong 
reduction of the Ncel temperature with doping; 

(e) a metal-insulator transition, driven by the tem- 
perature, takes place within the antiferromagnetic 
phase; at half filling and higher temperatures 
this transition coincides with the paramagnetic- 
antiferromagnetic one; 



(f) away from half-filling a metal-insulator transition 
driven by the doping is observed. This transition 
has the following properties: at zero temperature it 
is discontinuous and connects an antiferromagnetic 
metal and an antiferromagnetic insulator; at fi- 
nite temperature it connects the antiferromagnetic 
metal to an antiferromagnetic state of semiconduc- 
tor type. 

All these properties emerge from very strong correlations 
and are not usually found by approximations of mean- 
field type. They could be explained by analyzing the 
electronic density of states and the energy spectra. 

Finally, we want to point out that the thermodynam- 
ical study of the Hubbard model within the framework 
of the two-pole approximation by means of the COM is 
far from being completed. Besides the normal and the 
antiferromagnetic phase, a detailed analysis of the su- 
perconducting,^^ the ferromagnetic, the charge ordered 
and of other phases with more complex magnetic and 
charge ordering has to be completed or undertaken and 
these analysis should finally combine to give the com- 
plete phase diagram of the Hubbard model under this 
approach. 
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APPENDIX A: CALCULATION OF THE 
GREEN'S FUNCTIONS FOR THE t-t'-U MODEL 

On the magnetic lattice defined in Fig. 1 we 
can define^the translational Jnvariant Green's func- 
tions SJ^^(R, ,Rj,t),S^^{R^,Rj,t), 5*^^ (R^ , , t ) and 
S^^ (RijTLjjt) connecting two points Rj and R^ of the 
magnetic lattice. Their definition and the correspond- 
ing equations of motion arc illustrated in Fig. 21 for the 
translational invariant Green's function S^^. 
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FIG. 21. The definition of the translational invariant 
Green's function S^^. 

By means of the Fourier transform on the magnetic 
lattice 
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d^fc e'^<^'-^i^S^^{\i,uj) , X,Y€ {A, B} 



(1.1) 



the equations of motion for the translational invariant Green's functions take the form 

ujS^^{k,uj) = [/(") +/M] + ([^(i) +£(3)] +^(k)[£(5) +£(6)])5^A(k,^) + (5(k))*[£(2) +£W]5^^(k,a;) 
a;S^^(k,a;) = ([^(i) +^(3)] +^(k)[£(5) +e(6)])5AB(j^^^) ^ (5(k))*[e(2) +£(^)]5^^(k,a;) 

W5^^(k, = ([e(l) - £(3)] + ^(k) [e(5) _ ^(6)] \ gBAf^^^ ^) ^ ^^(2) _ ^(4)] ^^^(j^^ ^) 



a;5^^(k,a;) = [/(")-/«] + ([^d) -^O)] +/3(k)[£ 



r(5) _ ^(6) 



I 1 (k, o;) + 5(k) [£(2) _ £(4)] ^) 



(1.2) 



where we used the projections on the magnetic lattice 



5(k) = i[l 



i(4|5(k)|2 - 1) from Eq. (2.3). The expressions (3.17) 
and (2.1) are obtained from Eqs. (1.2) by lengthy but 
straightforward algebraic manipulations. 
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APPENDIX II: EXPLICIT EXPRESSIONS FOR 
THE GREEN'S FUNCTIONS OF THE 
ANTIFERROMAGNETIC EQUILIBRIUM STATE 

For X,Y e {A,B} the coefficients A^^ , B^^ , C^^ 
occurring in the expression (3.17) of the Green's func- 
tions take the explicit form: 



lAA 



j^BB 



(jBB 
j^BB 



O2 <- O2 
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^•BA ^ jiBB 
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(2.1) 



To abbreviate the notation we used the following defini- 
tions 



f+=£(i)+e(3) 

£t=e(^)+e^^) 
£:3+=£(5)+£(6) 



£1 

£- = £(2) _ £(4) 



£(5) 



£(3) 

(■ 

£(6) 



(2.2) 



Furthermore, we used the convention that cos (Q • Ri) 
takes positive values on the sublattice A, and wc notice 
that the projections on the nearest neighbors and on the 
next- and next-next-nearest neighbors according to the 
spherical approximation are given on the magnetic lattice 
by the expressions 



1 + e' 



5(k) = 
^(k) =i (e^^==^ 



j,o _|_ ^i{kxa+kya)'^ 
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i{k: 
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(2.3) 
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